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We are studying the effects of light dynamical quarks on the excitation energies of a flux tube between a static 
quark and antiquark. We report preliminary results of an analysis of the ground state potential and the Eg + 
and Il u potentials. We have measured these potentials on closely matched ensembles of gauge configurations, 
generated in the quenched approximation and with 2+1 flavors of Asqtad improved staggered quarks. 



1. INTRODUCTION 

Simulations with dynamical quarks have found 
that light quarks modify the heavy quark- 
antiquark potential in a number of ways Jl],||,[| ■ 
At large distances they decrease the string ten- 
sion in units of the Sommer r$ and r\ parame- 
ters (defined by r 2 F{r) = 1.65 and 1.00, respec- 
tively) and lead eventually to string-breaking. At 
shorter distances they modify the running of the 
coupling constant, deepening the Coulomb well 
and increasing the ratio ro/ri. In this work, we 
extend these studies to some of the potentials 
with excited flux tubes. Of particular interest to 
quarkonium spectroscopy are the H u excitations 
leading to exotic QQg hybrids j| . 

We report results of a study in which our 
sources and sinks are optimized to create and an- 
nihilate a flux-tube state. In the presence of dy- 
namical quarks, string breaking is expected. It is 
known that in the conventional channel, tran- 
sitions to the open two-meson channel are exceed- 
ingly weak, qualitatively consistent with the small 
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widths of quarkonium states above the heavy- 
light meson thresholds ||,||. Since at present we 
do not include the open two-meson channel we do 
not expect to observe string breaking here. 

2. MEASUREMENTS 

We have measured the heavy quark potential 
on an ensemble of 28 3 x 96 (a ~ 0.09 fm) gauge 
configurations generated in the presence of 2 + 1 
flavors of Asqtad dynamical quarks of varying 
masses and a one-loop Symanzik gauge action^. 
The strange quark mass is set approximately to 
its physical value. Here we compare results from 
our 358-configuration quenched ensemble with 
our 495-configuration dynamical quark ensemble 
for which (m w r ) 2 ~ 1.3. 

The configurations arc first smoothed using a 
single hypercubic (HYP) blocking pass 0, a tech- 
nique that improves significantly the signal-to- 
noisc ratio ||. The blocking procedure involves 
replacing all gauge links (timclikc as well as space- 
like) with an SU(3)-projected average over paths 
confined to adjacent hypercubes. Thus distor- 
tions in the result are local and expected to be 
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Figure 1. The ground state static quark po- 
tential for quenched (octagons) and 2+1 flavor 
(diamonds) QCD, in units of r . The solid lines 
are fits to the Coulomb plus constant plus linear 
form, fixing Vfit(ro) = 0. The lattice spacing is 
matched using r<j. The inset expands the area 
shown by the box. 



confined to distances smaller than about 2a g|. 
After HYP blocking the spacelike links are fur- 
ther smoothed via five cycles of APE smearing 
with SU(3) projection. 

On the thus smoothed lattices we measure the 
expectation value of the standard R x T Wilson 
loop on axis and along three different off-axis di- 
rections. These measurements yield the conven- 
tional ground state £+ and excited state E + ' po- 
tentials. For the n„ excited state, we measured 
the expectation value of a bent loop formed by 
replacing the source and sink flux tubes of length 
R by a superposition of large "staples" of sides 
(2a, R, 2a). For example one such loop replaces 
each on-axis spacelikc flux path (Rx) by paths of 
the form (2ay, Rx, —2ay) minus its reflection in 
the xz plane B. 
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Figure 2. The excited Y,' g potential for quenched 
(octagons) and 2+1 flavor (diamonds) QCD, in 
units of ro- The lattice spacing was matched us- 
ing r . 

For the standard Wilson loop we extracted the 
usual potential Vs g + and its excited state 
V^ g+ by doing a blocked, correlated, double- 
exponential fit to the Wilson loop data: 



W(R,T) = C Sg+ (R)e- y ^+W T 
+ C^ g+ {R)e- v ^ T . 



(1) 
(2) 



For the H u potential we did only a single- 
exponential fit. 

In all cases we use the same fit ranges for both 
quenched and dynamical lattices to reduce possi- 
ble systematic errors. 

3. RESULTS 

In Fig. |l| we compare the ground state poten- 
tial on the quenched ensemble and the 2+1 fla- 
vor ensemble. Both the distance scale and the 
potential arc plotted in units of ro, and a con- 
stant has been subtracted from the potential so 
that it is zero at ro. Since ro was determined 
from this potential, the fits are tangent at this 
point. Away from ro, the potentials have dif- 
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Figure 3. The n~ potential for quenched (oc- 
tagons) and 2+1 flavor (diamonds) QCD, in units 
of ro- The potentials are plotted relative to the 
zero determined in the fit to the ground state po- 
tentials. 



ferent shapes, namely, the Coulomb attraction is 
slightly stronger in the light quark ensemble and 
the string tension is slightly weaker in units of ro, 
confirming earlier findings Q . A softening of the 
Coulomb well is also evident. This is an expected 
consequence of HYP smoothing. 

Similarly in Fig. ^ we show light quark effects 
in the T,' g excitation potential. The potentials are 
plotted relative to the zero determined in the fit 
to the ground state potentials. From Fig. || we 
see that the excited state potential is slightly 
steeper than that of quenched QCD. 

In Fig. H we plot the H~ potential. This hybrid 
potential is weakly repulsive at short range, as 
would be expected from the Coulomb interaction 
in a color octet quark-antiquark system. This ef- 
fect is softened by HYP smoothing. At the level 
of our statistical errors there are no apparent dif- 
ferences at long range, but better statistics would 
certainly be of interest. 



4. CONCLUSIONS 

Our measurements at a = 0.09 fm confirm the 
shape changes in the ground-state potential, seen 
previously at a = 0.13 fm. In units of ro we 
find, further, that adding 2+1 flavors of dynam- 
ical quarks makes the E +/ excited state potential 
slightly steeper and the n~ slightly more repul- 
sive at short range. We find no clear evidence for 
a flattening of the potentials that would signal 
string breaking. 

Computations were performed at LANL, 
NERSC, NCSA, ORNL, PSC, SDSC, FNAL, and 
the CHPC (Utah). This work is supported by the 
U.S. NSF and DOE. 
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